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ABSTRACT 

We use multifrequency matched filters to estimate, in the WMAP 5-year data, the Sunyaev-Zel'dovich (SZ) fluxes of 893 ROSAT 
NORAS/REFLEX clusters spanning the luminosity range Lx.[o.i-2.4] kcv = 2x lO'" -3.5x lO"*^ erg s"' . The filters are spatially optimised 
by using the universal pressure profile recently obtained from combining XMM-Newton observations of the REXCESS sample and 
numerical simulations. Although the clusters are individually only marginally detected, we are able to firmly measure the SZ signal 
(> 10 o") when averaging the data in luminosity/mass bins. The comparison between the bin-averaged SZ signal versus luminosity 
and X-ray model predictions shows excellent agreement, implying that there is no deficit in SZ signal strength relative to expectations 
from the X-ray properties of clusters. Using the individual cluster SZ flux measurements, we directly constrain the F500 - i-x and 
K500 - M500 relations, where Ksoo is the Compton y-parameter integrated over a sphere of radius rjoo- The Ksoo-A^soo relation, derived 
for the first time in such a wide mass range, has a normalisation Kj^^ = [1.60 ± 0.19] x 10"^ arcmin" at Mjoo = 3 x lO''' /z"' Mq, in 
excellent agreement with the X-ray prediction of 1.54 x lO""* arcmin^, and a mass exponent of a = 1.79 ± 0.17, consistent with the 
self-similar expectation of 5 /3 . Constraints on the redshift exponent are weak due to the limited redshift range of the sample, although 
they are compatible with self-similar evolution. 

Key words. Cosmology: observations. Galaxies: cluster: general. Galaxies: clusters: intracluster medium. Cosmic background radi- 
ation. X-rays: galaxies: clusters 



1. Introduction parameter is given by the integral of the electron pressure along 

the line-of-sight in the direction h. 

Capability to observe the Sunyaev-Zel'dovich (SZ) effect has 

improved immensely in recent years. Dedicated instruments C kTg, 

y= -n.o-Tdl, (1) 
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now produce high resolution images of si ngle objects (e.g . 
iKitavama etal] 12004 iHalverson etaTI 12009- Nord et al. 2009) 



and moderately large samples of hig h-quaUty SZ measurements where ctt is the Thomson cross section. 

of previously-kno wn clusters (e.g., I Mroczkowski et alJ [2009t Most notably, the integrated SZ flux from a cluster directly 

iPlagge et al]|2009l) . In addition, large-scale surveys for clusters measures the total thermal energy of the gas. Expressing this flux 

using the SZ e ffect are underway, both from space with the in terms of the integrated Compton y-parameter Fsz- defined by 

Planck mission dValenziano et al. 1 12007: 'Lam arre et al. I l2003h j d£lAiy{h) = Fsz7y(-«) - we see that Fsz <^ J dndlrieT^ oc 

and from the ground with several dedicated telescopes, such j n.TedV. For this reason, we expect Fsz to closely correlate 

as flie South Pole Telescope dCarlstrom etalj|2009|) leading to ^j^h ^^^^1 cluster mass, M, and to provide a low-scatter mass 

the first discoveries of clusters solely through their SZ signal proxy 

dliniszewski et al. 2009). These results open the way for a bet- This expectatio n, bor ne out by b oth numerical simul ations 

ter understandmg of the 5Z - Mass relation and, ultimately, for (e.g.. lda Silva etZll^l^fotr^l2Mllj^ et al.ll2006h 

cosmological studies with large SZ cluster catalogues. indirectly from X-ray observations using Yx, the prod- 

The SZ effect probes the hot gas in the intracluster medium uct of the gas^ mass a nd mean temperature (Nagai et al. 200^ 

(ICM). Inverse Compton scattering of cosmic microwave back- lArnaud et alJ2007l : IVikhlinin et a l. 2009), strongly motivates the 

ground (CMB) photons by free electrons in the ICM creates a use of SZ cluster surveys as cosmological probes. Theory pre- 

unique spectral distortion dSunvaev & Zeldovich 1 [T970l[T97l diets the cluster abundance and its evolution - the mass func- 

seen as a frequency-dependent change in the CMB surface tion - in terms of M and the cosmological parameters. With 

brightness in the direction of galaxy clusters that can be writ- a good mass proxy, we can measure the mass function and its 

ten as A/v(n) - y{n) jyix), where jy is a universal function of evolution and hence constrain the cosmological model, includ- 

the dimensionless frequency x = hv/kTcmh- The Compton y- ing the properties of dark energy. In this context the relationship 
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between the integrated SZ flux and total mass, Fsz - M, is fun- 
damental as the required link between theory and observation. 
Unfortunately, despite its importance, we are only beginning to 
observationally cons train the relation dBonamente etal.ll2008t 
iMarrone et ani2009l) . 

Several a uthors have extracted the cluster SZ signal from 
WMAP data (iBennett et all I200I iHinshaw et ai]|2007l |2009|) . 
However, the latter are not ideal for SZ observations: the 
instrument having been designed to measure primary CMB 
anisotropies on scales larger than galaxy clusters, the spa- 
tial resolution and sensitivity of the sky maps render clus- 
ter detection difficult. Nevertheless, these authors extracted 
the cluster SZ signal by e ither cross-correlat i ng with the 
gener a l galaxy distribution ( Fosalba et al. 12003 : iMvers et al] 



l2004t iHernandez- Monteagudo et al.l |2004 l2006h or 'stack- 
ing' existing cl uster catalogues in t he optical or X-ray 



(Lieu et al. 2006; Afshordi et al. 2007; Atrio-Barande la et alj 
[2OO8; Bielbv & Shanks 2007; Diego & Partridge 2009) . These 
analyses indicate that an isothermal yS-model is not a good de- 
scription of the SZ profile, and some suggest that the SZ signal 
strength is lower th an exp ected from the X-ray properties of the 
clusters (iLieu et al.1l2006l:lBielbv & Shanksll2007l) . 

Recent in-depth X-ray studies of the ICM pressure profile 
demonstrate regularity in shape and simple scaling with clus- 
ter mass. Combining these observations with numerical simu- 
lations leads to a u niversal pressure profile (Nagai et al. 2007; 
lArnaud et al]|2009l) fliat is best fit by a modified NFW profile. 
The isothermal yS-model, on the other hand, does not provide an 
adequate fit. From this newly determined X-ray pressure pro- 
file, we can infer the expected SZ profile, y{r), and the Fsz - M 
relation at low redshift ( Arnaud et al. 200Sj). 

It is in light of this recent progress from X-ray obser- 
vations that we present a new analysis of the SZ effect in 
WMAP with the aim of constraining t he SZ scaling laws . 
We build a m ultifrequency matched filter (iHerranz et al. Il2002t 
iMelin et al.l2006i) based on the known spectral shape of the ther- 
mal SZ effect and th e shape of the universal pressure profile of 
[Arnaud et al.' ('2009'). This profile was derived from REXCESS 
(Bohringer et al. 2007), a sample expressly designed to measure 
the structural and scaling properties of the local X-ray cluster 
population by means of an unbiased, representative sampling in 
luminosity. Using the multifrequency matched filter, we search 
for the SZ effect in WMAP from a catalogue of 893 clusters de- 
tected by ROSAT, maximising the signal-to-noise by adapting 
the filter scale to the expected characteristic size of each clus- 
ter The size is estimated through the luminosity- mass relation 
derived from the REXCESS sample bv lPratt et all (12009). 

We then use our SZ measurements to directly determine the 
Ysz-Lx and Ys,z-M relations and compare to expectations based 
on the universal X -ray pressure profi l e. As comp ared to the pre- 
vious analyses of iBonamente et alTI (12008 ') and 'Marron e et al.1 
( 120091) . the large number of systems in our WMAP/ROS AT sam- 
ple allows us to constrain both the normalisation and slope of the 
Ysz - Lx and Fsz - M relations over a wider mass range and in 
the larger aperture of rsoo. In addition, the analysis is based on a 
more realistic pressure profile than in these analyses, which were 
based on an isothermal j0-model. Besides providing a direct con- 
straint on these relations, the good agreement with X-ray predic- 
tions implies that there is in fact no deficit in SZ signal strength 
relative to expectations from the X-ray properties of these clus- 
ters. 

The discussion proceeds as follows. We first present the 
WMAP 5-year data and the ROSAT cluster sample used, a 
combination of the REFLEX and NORAS catalogues. We then 



present the SZ model based on the X-ray-measured pressure 
profile (Sec.O. In Sec.|4] we discuss our SZ measurements, af- 
ter first describing how we extract the signal using the matched 
filter Section |5] details the error budget. We compare our mea- 
sured scaling laws to the X-ray predictions in Sect. |6] and Q and 
then conclude in Sec.|8] Finally, we collect useful SZ definitions 
and unit conversions in the Appendices. 

Throughout this paper, we use the WMAP5-only cosmolog- 
ical parameters set as our 'fiducial cosmology', i.e. h - 0.719, 
Qm = 0.26, Qa = 0.74, where h is the Hubble parameter at red- 
shift z = in units of 100 km/s/Mpc. We note /170 = /i/0.7 
and E{z) is the Hubble parameter at redshift z normalised to 
its present value. M500 is defined as the mass within the ra- 
dius r5oo at which the mean mass density is 500 times the 
critical density, Pcruiz), of the universe at the cluster redshift: 
M500 = ^TTpcririz) 500 rl^. 



2. The WMAP-5yr data and the NORAS/REFLEX 
cluster sample 

2.1. The WMAP-Syrdata 

We work with the WMAP full resolution coadded five year 
sky temperature maps at each frequency channel (downloaded 
from the LAMBDA archival). There are five full sky maps cor- 
responding to frequencies 23, 33, 41, 61, 94 GHz (bands K, 
Ka, Q, V, W respectively). The corresponding beam full widths 
at half maximum are approximately 52.8, 39.6, 30.6 21.0 and 
13.2 arcmins. The maps are orig inally at HEALPi)(Q resolution 
nside=512 (pixel= 6.87 arcmin). Although this is reasonably ad- 
equate to sample WMAP data, it is not adapted to the multi- 
frequency matched filter algorithm we use to extract the cluster 
fluxes. We oversample the original data, to obtain nside=2048 
maps, by zero-padding in harmonic space. In detail, this is per- 
formed by computing the harmonic transform of the original 
maps, and then performing the back transform towards a map 
with nside=2048, with a maximum value of { of /"^ax - 750, 850, 
1 100, 1500, 2000 for the K, Ka, Q, V, W bands respectively. The 
upgraded maps are smooth and do not show pixel edges as we 
would have obtained using the HEALPix upgrading software, 
based on the tree structure of the HEALPix pixelisation scheme. 
This smooth upgrading scheme is important as the high spatial 
frequency content induced by pixel edges would have been am- 
plified through the multifrequency matched filters implemented 
in harmonic space. 

In practice, the multifrequency matched filters are imple- 
mented locally on small, flat patches (gnomonic projection on 
tangential maps), which permits adaptation of the filter to the lo- 
cal conditions of noise and foreground contamination. We divide 
the sphere into 504 square tangential overlapping patches (100 
deg^ each, pixel=1.72 arcmin). All of the following analysis is 
done on these sky patches. 

The implementation of the matched filter requires knowl- 
edge of the WMAP beams. In this work, we assume symmet- 
ric beams, for which the transfer function b( is computed, in 
each frequency channel, from the noise-weighted average of the 
transfer functions of individual differential assemblies (a similar 
approach was used in Delab rouille et al...2009.) . 



http://lambda.gsfc.nasa.gov/ 
http://healpix.jpl.nasa.gov 
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Fig. 1. Inferred masses for the 893 NORAS/REFLEX clusters as 
a function of redshift. The cluster sample is flux limited. The 
right vertical axis gives the corresponding X-ray luminosities 
scaled by Eizy^^. The dashed blue lines del ineate th e mas s 
range over which the Lsoo-Msoo relation from jPratt et al] 120091) 
was derived. 



2.2. The NORAS/REFLEX cluster sample and derived X-ray 
properties 

We construct our cluster sample from the largest published X- 
ray catalogues: NORA S (^ohringer et al. 2000) and REFLEX 
dBohiringer et al.ll2004l) . both constructed from the ROSAT All- 
Sky S urvey. We merge the cluster lists given in Tabl es 1 , 6 and 8 
from Bohringer et al.l (2000) and Table 6 from Bo hringer et al.l 
(12004 ) and s nee the luminosities of the NORAS clusters are 
given in a standard cold dark matter (SCDM) cosmology {h - 
0.5, Qm = 1), we converted them to the WMAP5 cosmology. 
We also convert the luminosities of REFLEX clusters from the 
basic ACDM cosmology (h = 0.7, Qm = 0.3, Q.a = 0.7) to the 
more precise WMAP5 cosmology. Removing clusters appearing 
in both catalogues leaves 921 objects, of which 893 have mea- 
sured redshifts. We use these 893 clusters in the analysis detailed 
in the next Section. 

The NORAS/REFLEX luminosities Lx, measured in the soft 
[0.1-2.4] keV energy band, are given within various apertures 
depending on the cluster We convert the luminosities Lx to 
L500, the luminosities within rsoo, using an iterative scheme. 
This scheme is based on the mean electron density profile of 
the REXCESS cluster sample (Croston et al. 2008), which al- 
lows conversion of the luminosity be tween various ap ertures, 
and the REXCESS L500-M500 relation (iPratt et al.ll2009l) . which 



implicitly relates rsoo and Lsoo- The procedure thus simultane- 
ously yields an estimate of the cluster mass, M500, and the corre- 
sponding angular extent O^m = r5()o/Dang(z), where Dangiz) is the 
angular distance at redshift z. In the following we consider val- 
ues derived from relations both corrected and uncorrected for 
Malmquist bias. The relations are described by the following 
power law model^ 



E{z) 



-7/3 



-500 



C 



M 



M5, 



:00 



3 X 1014/!-(5M< 



(2) 



^ Since we consider a standard self-si milar model, we us ed the power 
law relations given in A ppendix B of lAmaud et alj ( l2009i) . They are 
derived as in iPratt et alj (2009) with the same luminosity data but for 
masses derived from a standard slope Msoo-ix relation. 



Table 1. Values for th e parameters of the Lx -M relation de rived 
from REXCESS data (iPratt et aI.ll2009HArnaud et al.ll2009l) 



Corrected for MB 


^°8(l0-l4;,-„^[ergs-"j 




Cr logL-logM 


no 


0.295 


1.50 


0.183 


yes 


0.215 


1.61 


0.199 



where the normalisation Cm, the exponent um and the dispersion 
(nearly constant with mass) are given in Table [1] The L500-M500 
relation was derived in the mass range [lO'^^-lO'^] M©. These 
limits are shown in Fig. [1] Note that we assume the relation is 
valid for lower masses. 

The final catalogue of 893 objects contains the position of 
the clusters (longitude and latitude), the measured redshift z, the 
derived X-ray luminosity L500, the mass M500 and the angular 
extent 6500. The clusters uniformly cover the celestial sphere at 
Galactic latitudes above \b\ > 20deg. Their luminosities L500 
range from 0.002 to 35.0 10**erg/s, and their redshifts from 
0.003 to 0.460. Figure[T]shows the masses M500 as a function of 
redshift z for the cluster sample (red crosses). The correspond- 
ing corrected luminosities L500 can be read on the right axis. The 
typical luminosity correction from measured Lx to L500 is about 
10%. The progressive mass cut-off" with redshift (only the most 
massive clusters are present at high z) reflects the flux limited 
nature of the sample. 

3. The cluster SZ model 

In this Section we describe the cluster SZ model, based on X-ray 
observations of the REXCESS sample combined with numer- 
ical simulations, as presented in Arnaud et al. (2009). We use 
the standard self-similar model presented in their Appendix B. 
Given a cluster mass M500 and redshift z, the model predicts the 
electronic pressure profile. This gives both the SZ profile shape 
and F500, the SZ flux integrated in a sphere of radius r^oQ. 



3.1. Cluster shape 

The dimensionless u niversal pressure p rofile is taken from 
Eq. Bl and Eq. B2 of lArnaud etaP (l2009l) : 



P(r) ^ Pq 

P500 xT'(l +x«)<^-T')/«' 



(3) 



wherex = r/r, withr, - rsoo/csoo andcsoo - 1.156, a - 1.0620, 
/3 = 5.4807, y = 0.3292 and with P500 defined in Eq.|4]below. 

This profile shape is used to optimise the SZ signal detec- 
tion. As described below in Sect.H] we extract the Ysz flux from 
WMAP data for each ROSAT system fixing C500, a, /3, j to the 
above values, but leaving the normalisation free. 

3.2. Normalisation 

The model allows us to compute the physical pressure profile as 
a function of mass and z, thus the Ysz-M^qo relation by integra- 
tion of P{r) to rgoo. For the shape parameters given above, the 
normalisation parameter Pq = 8.130 /Zj^'^ = 7.810 and the self- 
similar definition of f goo (lArnaud et al.ll2009l Eq. 5 and Eq. B2), 



P500 = 1.65xl0-'£(z) 



8/3 



M- 



500 



3 X lOl^/z-jjMoy 



2/3 



/i70 keV cm , 



(4) 
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Fig. 2. Left: Estimated SZ flux from a cylinder of aperture radius 5 x r^oo (Y^].^qq) as a function of the X-ray luminosity in an aperture 
of rjoo (Lsoo), for the 893 NORAS/REFLEX clusters. The individual clusters are barely detected. The bars give the total 1 cr error. 
Right: Red diamonds are the weighted average signal in 4 logarithmically-spaced luminosity bins. The two high luminosity bins 
exhibit significant SZ cluster flux. Note that we have divided the vertical scale by 30 between Fig. left and right. The thick and 
thin bars give the 1 cr statistical and total errors, respectively. Green triangles (shifted up by 20% with respect to diamonds for 
clarity) show the result of the same analysis when the fluxes of the clusters are estimated at random positions instead of true cluster 
positions. 



0.010 0.100 1.000 10.000 
L500 (0.1-2.4 keV) [10" erg/s] 



100.000 



one obtains: 



F500 [arcmin ] 



^500 



M500 



3 X 1014/z-lMe 



5/3 



E(z) 



2/3 



is covered by the five WMAP maps M,(jr) at frequencies v,=23, 
2 33, 41, 61, 94 GHz (/ = 1, ...,5). We arrange the survey maps 

igto a column vector M(x) whose component is the map at 

500 Mpc / ' frequency v,. The maps contain the cluster SZ signal plus noise: 



where F^gy - 1.54 x 10 (0T19) arcmin^. Equivalently, one 
can write: 



7500 [Mpc^] = Fg*oo 



M500 



5/3 



Eiz) 



2/3 



(6) 



where Y*^^ = 3.27 x 10"^ ((PTig) ^' Mpc^. Details of unit con- 
versions are given in Appendix iBl The mass dependence (M^^g) 

and the redshift dependence (Eizf'^^) of the relation are self- 
similar by construction. This model is used to predict the Y500 
value for each cluster These predictions are compared to the 
WMAP-measured values in Figs.[3]|ll|5]and|6] 



4. Extraction of the SZ flux 

4.1. Multifrequency Matched Filters 

We use multifrequency matched filters to estimate cluster fluxes 
from the WMAP frequency maps. By incorporating prior knowl- 
edge of the cluster signal, i.e., its spatial and spectral character- 
istics, the method maximally enhances the signal-to-noise of a 
SZ cluster source by optimally filtering the data. The universal 
profile shape described in Sec. [3] is assumed, and we evaluate 
the effects of uncertainty in this profile as outlined in Sec. |5] 
where we discuss our overall error budget. We fix the position 
and the characteristic radius 6^ of each cluster and estimate only 
its flux. The position is taken from the NORAS/REFLEX cat- 
alogue and 0s - S50o/c5()o with 0500 derived from X-ray data 
as described in Sec. 12.21 Below, we recall the main features of 
the multifrequency matched filters. More details can be found 
in|Herranz et al. (2002) or Melin et al. (2006). 

Consider a cluster with known radius 6^ and unknown central 
y-value yo positioned at a known point jCo on the sky. The region 



M(x) = yojvTgXx - Xo) + N(x) 



(7) 



where is the noise vector (whose components are noise maps 
at the different observation frequencies) and jV is a vector with 
components given by the SZ spectral function jy evaluated at 
each frequency. Noise in this context refers to both instrumental 
noise as well as all signals other than the cluster thermal SZ ef- 
fect; it thus also comprises astrophysical foregrounds, for exam- 
ple, the primary CMB anisotropy, diffuse Galactic emission and 
extragalactic point sources. Tg^(x - Xo) is the SZ template, tak- 
ing into account the WMAP beam, at projected distance (x - Xo) 
from the cluster centre, normalised to a central value of unity 
before convolution. It is computed by integrating along the line- 
of-sight and normalising the universal pressure profile (Eq.[3]l. 
The profile is truncated at 5 x rsoo (i.e. beyond the virial radius) 
so that what is actually measured is the flux within a cylinder of 
aperture radius 5 x rjoo- 

X-ray observations are typically well-constrained out to rsoo. 
Our decision to integrate out to 5 x rsoo is motivated by the fact 
that for the majority of clusters the radius rsoo is of order the 
Healpix pixel size (nside=512, pixel=6.87 arcmin). Integrating 
only out to rsoo would have required taking into account that 
only a fraction of the flux of some pixels contributes to the true 
SZ flux in a cylinder of aperture radius rsoo- We thus obtain the 
total SZ flux of each cluster by integrating out to 5 x rsoo, and 
then convert this to the value in a sphere of radius rsoo for direct 
comparison with the X-ray prediction. 

The multifrequency matched filters ^g^ix) return a minimum 
variance unbiased estimate, fo, of yo when centered on the clus- 
ter: 



fo = j d\x^gJ(x-Xo)-M{x) 



(8) 



where superscript t indicates a transpose (with complex conju- 
gation when necessary). This is just a linear combination of the 
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maps, each convolved with its frequency-specific filter 
The result expressed in Fourier space is: 



where 

FeSk) = jvTeJk) 



dhFgJ(k)-p-' -FeXk) 



-1/2 



(9) 

(10) 
(11) 



with P(k) being the noise power spectrum, a matrix in fre- 
quency space with components Pjj defined by {Ni{k)N*{k')}N - 
Pij{k)6(k - k'). The quantity erg gives the total noise variance 
through the filter, corresponding to the statistical errors quoted 
in this paper The other uncertainties are estimated separately as 
described in Sec. 15. II The noise power spectrum P(k) is directly 
estimated from the maps: since the SZ signal is subdominant 
at each frequency, we assume N(x) ~ M{x) to do this calcula- 
tion. We undertake the Fourier transform of the maps and aver- 
age their cross-spectra in annuli with width Al = 180. 

4.2. Measurements of the SZ flux 

The derived total WMAP flux from a cylinder of aperture radius 

Sxrsoo (J's^soo) for *e 893 individual NORAS/REFLEX clusters 
is shown as a function of the measured X-ray luminosity L500 in 
the left-hand panel of Fig. |2l The clusters are barely detected 
individually. The average signal-to-noise ratio (S/N) of the total 
population is 0.26 and only 29 clusters are detected atS/N > 2, 
the highest detection being at 4.2. However, one can distinguish 
the deviation towards positive flux at the very high luminosity 
end. 

In the right-hand panel of Fig. |2] we average the 893 mea- 
surements in four logarithmically-spaced luminosity bins (red 
diamonds plotted at bin center). The number of clusters are 7, 
150, 657, 79 from the lowest to the highest luminosity bin. Here 
and in the following, the bin average is defined as the weighted 
mean of the SZ flux in the bin (weight of l/cTg ). The thick er- 
ror bars correspond to the statistical uncertainties on the WMAP 
data only, while the thin bar gives the total errors as discussed 
in Sec. 15. II The SZ signal is clearly detected in the two highest 
luminosity bins (at 6.0 and 5.4 cr, respectively). As a demonstra- 
tive check, we have undertaken the analysis a second time using 
random cluster positions. The result is shown by the green trian- 
gles in Fig.|2]and is consistent with no SZ signal, as expected. 

In the following Sections, we study both the relation be- 
tween the SZ signal and the X-ray luminosity and that with the 
mass Msoo- We consider Ysm, the SZ flux from a sphere of ra- 
dius rgoo, converting the measured Y'^-'J^^^^ into Ysoo as described 
in Appendix |A] This allows a more direct comparison with the 
model derived from X-ray observations (Sec.|3]l. Before present- 
ing the results, we first discuss the overall error budget. 



5. Overall error budget 

5.1. Error due to dispersion in X-ray properties 

The error erg on Y500 given by the multifrequency matched fil- 
ter only includes the statistical SZ measurement error, due to the 
instrument (beam, noise) and to the astrophysical contaminants 
(primary CMB, Galaxy, point sources). However, we must also 
take into account: 1) uncertainties on the cluster mass estimation 



from the X-ray luminosities via the L500 - M500 relation, 2) un- 
certainties on the cluster profile parameters. These are sources of 
error on individual Fgoo estimates (actual parameters for each in- 
dividual cluster may deviate somewhat from the average cluster 
model). These deviations from the mean, however, induce addi- 
tional random uncertainties on statistical quantities derived from 
Fsoo, i.e. bin averaged Y^m values and the ^500-^500 scaling rela- 
tion parameters. Their impact on the isoo-Mioo relation, which 
depends directly on the M500 estimates, is also an additional ran- 
dom uncertainty. 

The uncertainty on M500 is dominated by the intrinsic disper- 
sion in the Lsoo-Msoo relation. Its effect is estimated by a Monte 
Carlo (MC) analysis of 100 realisations. We use the dispersion 
at z = as estimated bv lPratt et al.l (l2009l) . given in Table [T] For 
each realisation, we draw a random mass log M500 for each clus- 
ter from a Gaussian distribution with mean given by the Lsoo- 
M500 relation and standard deviation cTiagi^iogM/aM- We then 
redo the full analysis (up to the fitting of the Ysz scaling rela- 
tions) with the new values of M500 (thus 0s)- 

The second uncertainty is due to the observed dispersion 
in the cluster pr ofile shape, which depends on radius as shown 
in lArnaud et al.l l2009i (T\ogp ~ 0.10 beyond the core). Using 
new 100 MC realisations, we estimate this error by drawing a 
cluster profile in the log-log plane from a Gaussian distribu- 
tion with mean given by Eq. [3] and standard deviation depend- 
in g on the cluster radi us as shown in the lower panel of Fig. 2 
in lArnaud et al.l(l2009l) . 

The total error on Y^m and on the scaling law parameters 
is calculated from the quadratic sum of the standard deviation 
of both the above MC realisations and the error due to the SZ 
measurement uncertainty. 

5.2. Tlie Malmquist bias 

The NORAS/REFLEX sample is flux limited and is thus subject 
to the Malmquist bias (MB). This is a source of systematic er- 
ror Ideally we should use a ^500-^500 relation which takes into 
account the specific bias of the sample, i.e. computed from the 
true Lsoo-Msoo relation, with dispersion and bias according to 
each survey selection function. We have an estimate of the true, 
ie MB corrected, L500-M500 relation, from the published analy- 
sis of REXCESS data (Table 1). However, while the REFLEX 
selection function is known and available, this is not the case 
for the NORAS sample. This means that we cannot perform a 
fully consistent analysis. In order to estimate the impact of the 
Malmquist Bias we thus present, in the following, results for two 
cases. 

In the first case, we use the published Lsoo-Msoo relation 
derived directly from the REXCESS data, i.e. not corrected for 
the REXCESS MB (hereafter the REXCESS L50Q-M500 relation). 
Note that the REXCESS is a sub-sample of REFLEX. Using this 
relation should result in correct masses if the Malmquist bias 
for the NORAS/REFLEX sample is the same as that for the 
REXCESS. The YsoQ-MsQo relation derived in this case would 
also be correct and could be consistently compared with the X- 
ray predicted relation. We recall that this relation was derived 
from pressure and mass measurements that are not sensitive to 
the Malmquist bias. However Lsoo would remain uncorrected so 
that the Fsoo-Ljoo relation derived in this case should be viewed 
as a relation uncorrected for the Malmquist bias. In the second 
case, we use the MB corrected L500-M500 relation (hereafter the 
intrinsic Lsoo-Msoo relation). This reduces to assuming that the 
Malmquist bias is negligible for the NORAS/REFLEX sample. 
The comparison of the two analyses provides an estimate of the 
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direction and amplitude of the effect of the Malmquist bias on 
our results. The REXCESS Lsoo-A^soo relation is expected to be 
closer to the Lsoo-Msoo relation for the NORAS/REFLEX sam- 
ple than the intrinsic relation. The discussions and figures cor- 
respond to the results obtained when using the former, unless 
explicitly specified. 

The choice of the L500-M500 relation has an effect both on 
the estimated L500, Msoo and Y^oq values and on the expectation 
for the SZ signal from the NORAS/REFLEX clusters. However, 
for a cluster of given luminosity measured a given aperture, L500 
depends weakly on the exact value of rsoo due to the steep drop 
of X-ray emission with radius. As a result, and although Lsoo 
and M500 (or equivalently rsoo) are determined jointly in the it- 
erative procedure described in Sec. 12.21 changing the underly- 
ing Ljoo-Msoo relation mostly impacts the M500 estimate: L500 
is essentially unchanged (median difference of ~ 0.8%) and the 
difference in M500 simply reflects the difference between the re- 
lations at fixed luminosity. This has an impact on the measured 
F500 via the value of rsoo (the profile shape being fixed) but the 
effect is also small (< 1%). This is due to the rapidly converging 
nature of the Ysz flux (see Fig. 1 1 of Arnaud et al. 2009). On the 
other hand, all results that depend directly on Msoo, namely the 
derived Fsoo-A^soo relation or the model value for each cluster, 
that varies as M^^g (Eq. |5]), depend sensitively on the L500-M500 
relation. M500 derived from the intrinsic relation is higher, an ef- 
fe ct increasing wit h decreasing cluster luminosity (see Fig. B2 
of lPratt et al.l[2009h . 



5.3. Other possible sources of uncertainty 

The analysis presented in this paper has been performed on 
the entire NORAS/REFLEX cluster sample without removal of 
clusters hosting radio point sources. To investigate the impact 
of the p oint sources o n our result, we have cross-correlated the 
NVSS (iCondon et all 1 1998) and SUMMS (iMauch et al.ll2003l) 
catalogues with our cluster catalogue. We conservatively re- 
moved from the analysis all the clusters hosting a total radio 
flux greater than 1 Jy within 5 x r^m. This leaves 328 clusters in 
the catalogue, removing the measurements with large uncertain- 
ties visible in Figure |2] left. We then performed the full analysis 
on these 328 objects up to the fitting of the scaling laws, find- 
ing that the impact on the fitted values is marginal. For example, 
for the REXCESS case, the normalisation of the Fsoo-Msoo rela- 
tion decreases from 1.60 to 1.37 (1.6 statistical cr) and the slope 
changes from 1.79 to 1.64 (1 statistical cr). The statistical errors 
on these parameters decrease respectively from 0.14 to 0.30 and 
from 0.15 to 0.40 due to the smaller number of remaining clus- 
ters in the sample. 

The detection method does not take into account superposi- 
tion effects along the line of sight, a drawback that is inherent to 
any SZ observation. Thus we cannot fully rule out that our flux 
estimates are not partially contaminated by low mass systems 
surrounding the clusters of our sample. Numeri cal simulations 
give a possible estimate of the contamination: i Hallman et al.l 
(l2007h suggest that low-mass systems and unbound gas may 
contribute up to 16.3%^™^^ of the SZ signal. This would lower 
our estimated cluster fluxes by ~ 1.5cr. 



6. The Fsz"^5oo relation 

6.1. WMAP SZ measurements vs. X-ray model 

We first consider bin averaged data, focusing on the luminosity 
range Lsoo > lO**^ ergs/s where the SZ signal is significantly de- 
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Fig. 4. Estimated SZ flux F500 (in a sphere of radius rsoo) as a 
function of the mass Msoo averaged in 4 mass bins. Red dia- 
monds are the WMAP data. Blue stars correspond to the X-ray 
based model predictions and are shifted to higher masses by 20% 
for clarity. The model is in very good agreement with the data. 

tected (Fig. |2] right). The left panel of Fig. [3] shows Fsoo from a 
sphere of radius rsoo as a function of L500, averaging the data in 
six equally-spaced logarithmic bins in X-ray luminosity. Both 
quantities are scaled according to their expected redshift depen- 
dence. The results are presented for the analyses based on the 
REXCESS (red diamonds) and intrinsic (red crosses) Lsoo-Msoo 
relations. For the reasons discussed in Sec. 15.21 the derived data 
points do not differ significantly between the two analyses (Fig.|3] 
left), confirming that the measured Fgoo-isoo relation is insensi- 
tive to the finer details of the underlying Lsoo-M^m relation. 

We also apply the same averaging procedure to the model 
F500 values derived for each cluster in Sec. [3] The expected val- 
ues for the same luminosity bins are plotted as stars in the left- 
hand hand panel of Fig. [3] The Fsoo-isoo relation expected from 
the combination of the Fsoo-A^soo (Eq- HJ and L500-M500 (Eq. |2]) 
relations is superimposed to guide the eye. The right-hand panel 
of Fig. [3] shows the ratio between the measured data points and 
those expected from the model. As discussed in Sec. 15.21 the 
model values depend on the assumed L500-M500 relation. The 
difference is maximum in the lowest luminosity bin where the 
intrinsic relation yields ~ 40% higher value than the REXCESS 
relation (Fig. [3] left panel). The SZ model prediction and the 
data are in good agreement, but the agreement is better when 
the REXCESS L5oo-A^50o is used in the analysis (Fig. [3] right 
panel). This is expected if indeed the agreement is real and the 
effective Malmquist bias for the NORAS/REFLEX sample is not 
negligible and is similar to that of the REXCESS. 

6.2. 7500-^500 relation fit 

Working now with the individual flux measurements, F500, and 
L500 values, we fit an Fsoo-i'soo relation of the form: 

_ , / £(z)-^/^L5oo f / Da„,(z) 

using the statistical error on F500 given by the multifrequency 
matched filter The total error is estimated by Monte Carlo (see 
Sec. 15.11 ) but is dominated by the statistical error. The results 
are presented in Table|2] As already stated in Sec. 16. II the fitted 
values depend only weakly on the choice of Lsoo-Msoo relation. 
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Fig. 3. Left: Bin averaged SZ flux from a sphere of radius rsoo (Fsoo) as a function of X-ray luminosity in a aperture of rsoo (isoo)- 
The WMAP data (red diamonds and crosses), the SZ cluster signal expected from the X-ray based model (blue stars) and the 
combination of the Ysm-Msm and L500-M500 relations (dash and dotted dashed lines) are given for two analyses, using respectively 
the intrinsic L500-M500 and the REXCESS L^oo-Msa) relations. As expected, the data points do not change significantly from one 
case to the other showing that the Fjoo-isoo relation is rather insensitive to the finer details of the underlying Lsoo-Msoo relation. 
Right: Ratio of data points to model for the two analysis. The points for the analysis undertaken with the intrinsic Lsoo-Msoo are 
shifted to lower luminosities by 20% for clarity. 

Table 2. Fitted parameters for the observed Ysz-L^qq relation. The X-ray based model gives Y^^ = 0.89|1.07 x 

10"^ (/l/0.719)"^^^arcmin^ 

Oy - 1.1 1|1.04 and /By = 2/3 for the REXCESS and intrinsic L500-M500 relation, respectively. 



L500 - M500 K*^n [10-^ (h/0Jl9)-' arcmini 






REXCESS 0.92 ± 0.08 stat [+0.10 tot] 


1.11 (fixed) 


2/3 (fixed) 


0.88 +0.10 Stat [±0.12 tot] 


1.19 ±0.10 stat [±0.10 tot] 


2/3 (fixed) 


0.90 ±0.13 Stat [±0.16 tot] 


1.11 (fixed) 


1.05 ±2.18stat [±2.25 tot] 


Intrinsic 0.95 ± 0.09 stat [±0. 1 1 tot] 


1.04 (fixed) 


2/3 (fixed) 


0.89 ±0.10 Stat [±0.12 tot] 


1.19 ±0.10 stat [±0.10 tot] 


2/3 (fixed) 


0.89 ±0.13 stat [±0.16 tot] 


1.04 (fixed) 


2.06 ±2. 14 stat [±2.21 tot] 



7. The FSZ-M500 relation and its evolution 

In this Section, we study the mass and redshift dependence 
of the SZ signal and check it against the X-ray based model. 
Furthermore, we fit the isoo-Mgoo relation and compare it with 
the X-ray predictions. 

7. 1 . Mass dependence and redshift evolution 

Figure |4] shows the bin averaged SZ flux measurement as a func- 
tion of mass compared to the X-ray based model prediction. As 
expected, the SZ cluster flux increases as a function of mass and 
is compatible with the model. In order to study the behaviour of 
the SZ flux with redshift, we subdivide each of the four mass bins 
into three redshift bins coiTesponding to the following ranges: 
z < 0.08, 0.08 < z < 0.18, z > 0.18. The result is shown 
in the left panel of Fig. |5] In a given mass bin the SZ flux de- 
creases with redshift, tracing the £>a„g(z)"^ dependence of the 
flux. In particular, in the highest mass bin (lO'^M©), the SZ flux 
decreases from 0.007 to 0.001 arcmin^ while the redshift varies 
fromz < 0.08 toz > 0.18. The mass and the redshift dependence 
are in good agreement with the model (stars) described in Sec. [3] 
Since the DangizY^ dependence is the dominant effect in the 
redshift evolution, we multiply Fgoo by Dangiz)^ and divide it 



by the self-similar mass dependence M^qq- expected self- 
similar behaviour of the new quantity J'5oo£'flnj(z) /.W500 a 
function of redshift is E(z)^^^ (see Eq. |5]l. The right panel of 
Fig. |5] shows YsQoDangiz)^ /M^l^ as a function of redshift for the 
three redshift bins z < 0.08, 0.08 < z < 0.18, z > 0.18. The 
points have been centered at the average value of the cluster red- 
shifts in each bin. The model is displayed as blue stars. Since 
the model has a self-similar redshift dependence and Eiz)^^^ in- 
creases only by a factor of 5% over the studied redshift range, 
the model stays nearly constant. The blue dotted line is plotted 
through the model and varies as E{z)^^^. The data points are in 
good agreement with the model, but clearly, the redshift lever- 
age of the sample is insufficient to put strong constraints on the 
evolution of the scaling laws. 



We now focus on the mass dependence of the relation. We 
scale the SZ flux with the expected redshift dependence and plot 
it as a function of mass. The result is shown in Fig.|6]for the high 
mass end. The figure shows a very good agreement between the 
data points and the model, which is confirmed by fitting the rela- 
tion to the individual SZ flux measurements (see next Section). 
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Fig. 5. Evolution of the ^500-^^500 relation. Left: The WMAP data from Fig. |4]are divided into three redshift bins: z<0.08 (blue 
diamonds), 0.08<z<0.18 (green crosses), z>0.18 (red triangles). We observe the expected trend: at fixed mass, F500 decreases with 



redshift. This redshift dependence is mainly due to the angular distance (Fgoo Da„g{z) ^). The stars give the prediction of the 



model. Right: We divide F500 by M^^^Da„g(z) and plot it as a function of z to search for evidence of evolution in the Fjoo-A^soo 
relation. The thick bars give the 1 cr statistical errors from WMAP data. The thin bars give the total 1 sigma errors. 
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Table 3. Fitted parameters for the observed Ysz-M^qq relation. The X-ray based model gives Y'^^ 
10-^ (/i/0.719)"^^^arcmin^ ay = 5/3 andySy = 2/3. 



1.54 X 



L500-M500 relation Y*,„ [10"^ (A/0.719)"- arcmin^] 


ay 


/3y 


REXCESS 1.60 ±0.14 Stat [+0.1 9 tot] 


5/3 (fixed) 


2/3 (fixed) 


1.60 ±0.15 Stat [±0.19 tot] 


1.79 ±0.15 Stat [±0.17 tot] 


2/3 (fixed) 


1.57 ±0.23 Stat [±0.29 tot] 


5/3 (fixed) 


1.05 ± 2.18 Stat [±2.52 tot] 


intrinsic 1.37 ± 0.12 stat [±0.17 tot] 


5/3 (fixed) 


2/3 (fixed) 


1.36 ±0.13 Stat [±0.17 tot] 


1.91 ±0.1 6 Stat [±0.1 8 tot] 


2/3 (fixed) 


1.28 ±0.19 Stat [±0.24 tot] 


5/3 (fixed) 


2.06 ±2. 14 Stat [±248 tot] 



7.2. Ysoo-Msoo relation fit 

Using the individual F500 measurements and M500 estimated 
from the X-ray luminosity, we fit a relation of the form: 



Y500 



^500 



500 



3 X 



E(zf' 



Dgngiz) 

500 Mpc 



(13) 



The results are presented in Table [3] for the analysis undertaken 
using the REXCESS and that using the intrinsic Lsoo-Msoo re- 
lation. The p ivot mass 3 x lO'^'/z^'M©, close to that used by 
lAmaud et alJ (2009), is slightly larger than the average mass of 
the sample (2.8|2.5 x IO'^^Mq for the REXCESS — intrinsic L500- 
M500 relation, respectively). We use a non-linear least-squares 
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fit built on a gradient-expansion algorithm (IDL curvefit func- 
tion). In the fitting procedure, only the statistical errors given by 
the matched multifilter (crY^gJ are taken into account. The total 
errors on the final fitted parameters, taking into account uncer- 
tainties in X-ray properties, are estimated by Monte Carlo as 
described in Sec.|5] 

We first discuss the results obtained using the R EXCESS 
Lsoo-Msoo relation, which is expected to be close to the opti- 
mal case (see discussion in Sec. 15.2b . First, we keep the mass 
and redshift dependence fixed to the self-similar expectation 
(ay - 5/3, ySy = 2/3) and we fit only the normalisation. We ob- 
tain Y^QQ - 1.60 X lO"-' (/i/0.719)"^ arcmin^, in agreement with 

the X-ray prediction F'^^ = 1.54 x 10"^ (/j/0.719)"^^^ arcmin^ 
(at QAcr). Then, we relax the constraint on ay and fit the normal- 
isation and mass dependence at the same time. We obtain a value 
for ay = 1 .79, slightly greater than the self-similar expectation 
(5/3) by 0.8 cr. To study the redshift dependence of the effect, we 
fix the mass dependence to ay =5/3 and fit Fgy,, and /Sy at the 
same time. We obtain a somewhat stronger evolution fiy = 1.05 
than the self-similar expectation (2/3). The difference, however, 
is not significant (0.2cr). As already mentioned above (see also 
Fig. |5] right), the redshift leverage is too small to get interesting 
constraints on /3y. 

As cluster mass estimates depend on the assumption of the 
underlying L500-M500 relation, so does the derived Ysm-M^m 
relation as well. However, the effect is small. The normalisa- 
tion is shifted from (1.60 + 0.14 stat [+0.19 tot]) lO^^arcmin^ 
to (1.37 + 0.12 Stat [+0.17 tot]) 10"^ arcmin^ when using the in- 
trinsic Lsoo-A^soo relation. The difference is less than two statis- 
tical sigmas, and for the mass exponent, it is less than one. 



8. Discussion and conclusions 

In this paper we have investigated the SZ effect and its scal- 
ing with mass and X-ray luminosity using WMAP 5-year data 
of the largest published X-ray-selected cluster catalogue to 
date, derived from the combined NORAS and REFLEX sam- 
ples. Cluster SZ flux estimates were made using an optimised 
multifrequency matched filter. Filter optimisation was achieved 
through priors on the pressure distribution (i.e., cluster shape) 
and the integration aperture (i.e., cluster size). The pressure dis- 
tribution is assumed to follow the universal pressure profile of 
[Arnaud et al. ( 2009), derived from X-ray observations of the rep- 
resentative local REXCESS sample. This profile is the most re- 
alistic available for the general population at this time, and has 
been shown to be in good a greement with rec ent high-quality SZ 
observations from SPT (IPla^^e et alJl200^ . Furthermore, our 
analysis takes into account the dispersion in the pressure distri- 
bution. The integration aperture is estimated from the L500-M500 
relation of the same REXCESS sample. We emphasise that these 
two priors determine only the input spatial distribution of the 
SZ flux for use by the multifrequency matched filters; the priors 
give no information on the amplitude of the measurement. As 
the analysis uses minimal X-ray data input, the measured and 
X-ray predicted SZ fluxes are essentially independent. 

We studied the Ysz - Lx relation using both bin averaged 
analyses and individual flux measurements. The fits using indi- 
vidual flux measurements give quantitative results for calibrat- 
ing the scaling laws. The bin averaged analyses allow a direct 
quantitative check of SZ flux measurements versus X-ray model 
predictions based o n the universal pressure profile derived by 
lArnaud et"al] (l2009h from REXCESS. An excellent agreement is 
found. 



Using WMAP 3-ye ar data, both iLieu et al. I (l2006t) and 
iBielbv & Shanks! (12007!) found that the SZ signal strength is 
lower than predicted given expectations from the X-ray proper- 
ties of their clusters, concluding that that there is some missing 
hot gas in the intra-cluster medium. The excellent agreement be- 
tween the SZ and X-ray properties of the clusters in our sample 
shows that there is in fact no deficit in SZ signal strength relative 
to expectations from X-ray observations. Due to the large size 
and homogeneous nature of our sample, and the internal con- 
sistency of our baseline cluster model, we believe our results 
to be robust in this respect. We note that there is some con- 
fusion in the literature regarding the phrase 'missing baryon s' . 
The 'missing baryons' mentioned by !Afshordi et al.l 12007!) in 
the WMAP 3-year data are missing with respect to the univer- 
sal baryon fraction, but not with respect to the expectations from 
X-ray measurements. Afshordi et al. (2007) actually found good 
agreement between the strength of the SZ signal and the X-ray 
properties of their cluster sample, a conclusion that agrees with 
our results. This good convergence between SZ direct measure- 
ments and X-ray data is an encouraging step forward for the 
prediction and interpretation of SZ surveys. 

Using L500 as a mass proxy, we also calibrated the Fsoo-Mjoo 
relation, finding a normalisation in excellent agreement with X- 
ray predictions based on the universal pressure profile, and a 
slope consistent with self-similar expectations. However, there 
is some indication that the slope may be steeper, as also indi- 
cated from the REXCESS a nalysis when using the best fitting 
empirical Msoo-ix relation dArnaud et aljr2009 ). Msoo depends 
on the assumed L500 - M500 relation, making the derived F500- 
M500 relation sensitive to Malmquist bias which we cannot fully 
account for in our analysis. However, we have shown that the 
effect of Malmquist bias on the present results is less than 2cr 
(statistical). 

Regarding evolution, we have shown observationally that the 
SZ flux is indeed sensitive to the angular size of the cluster 
through the diameter distance effect. For a given mass, a low 
redshift cluster has a bigger integrated SZ flux than a similar 
system at high redshift, and the redshift dependence of the in- 
tegrated SZ flux is dominated by the angular diameter distance 
(oc D^„j(z) ). However, the redshift leverage of the present clus- 
ter sample is too small to put strong contraints on the evolution 
of the YsoQ-L^oo and Fsoo-A^soo relations. We have nevertheless 
checked that the observed evolution is indeed compatible with 
the self-similar prediction. 

In this analysis, we have compensated for the poor sensitivity 
and resolution of the WMAP experiment (regarding SZ science) 
with the large number of known ROSAT clusters, leading to self- 
consistent and robust results. We expect further progress using 
upcoming Planck all-sky data. While Planck will offer the possi- 
bility of detecting the clusters used in this analysis to higher pre- 
cision, thus significantly reducing the uncertainty on individual 
measurements, the question of evolution will not be answered 
with the present RASS sample due to its limited redshift range. 
A complementary approach will thus be to obtain new high sen- 
sitivity SZ observations of a smaller sample. The sample must 
be representative, cover a wide mass range, and extend to higher 
z (e.g., XMM-Newton follow-up of samples drawn from Planck 
and ground based SZ surveys). This would deliver efficient con- 
straints not only on the normalisation and slope of the Ysz - Lx 
and Ysz - M relations, but also their evolution, opening the way 
for the use of SZ surveys for precision cosmology. 
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Appendix A: SZ flux definitions 

In this Appendix, we give the definitions of SZ fluxes we used. 
Table lA.ll gives the equivalence between them. In this paper, we 
mainly use F500 as the definition of the SZ flux. This flux is the 
integrated SZ flux from a sphere of radius rsoo- It can be related 
to Y„r5Qo, the flux from a sphere of radius n x rsoo by integrating 
over the cluster profile: 
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drP(r)4nr^ 



Jg drf{r)^nr 
Jy drP{r)Anr 



(A.l) 



where P{r) is given by Eq. [3] The ratio Fnrsoo / i'soo is given in 
TableOfor « = 1, 2 , 3 , 5 , 10. 

In practice, an experiment does not directly measure Ysqq but 
the SZ signal of a cluster integrated along the line of sight and 
within an angular aperture. This corresponds to the Compton pa- 
rameter integrated over a cylindrical volume. In Sec. ID we esti- 
mate Y^^^^, the flux from a cylinder of aperture radius 5 x rsoo 
using the matched multifilter. Given the cluster profile, we can 
derive F„r=i()o from Y' 



cyl 
'«r500- 



JO Jr s: 



sm6<iir5QQ 



d6P{r)2nr^ 



P'" drP(r)4nr^ 



(A.2) 



«r500 



/y: 



nrSOO 



given in Table lA.ll for n = 



The ratio Y, 

1,2,3,5,10. In the paper, we calculate isoo from I500 = 
0.986/1.814x7; 



■cyl 

5/-500- 



Appendix B: SZ units conversion 

In this Appendix, we provide the numerical factor needed for 
the SZ flux units conversion and derive the relation between the 
recently introduced Yx parameter and the SZ flux Ysz- The latter 
will allow readers to easily convert between SZ fluxes given in 
this paper and those reported in other publications. 
Given the definition of SZ flux: 



' nrSOO 



-j 

Jn„ 



dQy 



(B.l) 



where Q„r500 is the solid angle covered by n x rsoo, and the fact 
that the Compton parameter y is unitless, the observational units 
for the SZ flux are those of a solid angle and usually given in 
arcmin^. 

The SZ flux can be also computed in units of Mpc^ and the 
conversion is given by 

y.z[Mpc^] . 60- (^J y.z[arcmin^] (^)' 



8.46 x 10"'^ Ysz[arcmm^] 



Dangjz) 

1 Mpc 



(B.2) 



Table A.l. Equivalence of SZ flux definitions 



n 


1 


2 


3 


5 


10 


YnriOolYsOO 


1 


1.505 


1.690 


1.814 


1.873 


^"'500/^5(10 


0.827 


0.930 


0.963 


0.986 


0.997 
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where Da„g(z) is the angular distance to the cluster. 

The X-ray analogue of the integrated SZ Comptonisation pa- 
rameter is Yx - Mgas,50()Tx whose natural units are keV, 
where Mgas^sm is the gas mass in rsoo and Tx is the spec- 
troscopic temperatur e excluding the central 0.15 rjoo region 
(iKravtsov et al.ll2006l) . To convert between Ysz and Yx, we first 
have 



f 

Jo 



Fsz[Mpc2]= I " t/rcrr ^n.WW (B.3) 



where crj is the Thomson cross section (in Mpc ), m^c the elec- 
tron mass (in keV), Ti.{r) the electronic temperature (in keV) and 
Heir) the electronic density. By assuming that the gas tempera- 
ture Tg(r) is equal to the electronic temperature Teif) and writ- 
ing the gas density as Pgif) - iiempii ,,{/), where nip is the proton 
mass and //^ = 1.14 the mean molecular weight per free electron, 
one obtains: 



(T 1 /'^''soo 

Fsz[Mpc2] ^ —^-—\ drpg(r)Tg(r)Anr^ 



1 r 

l^emp Jo 

= Cxsz A^ga.5,500 Tmw - A Cxsz (B-4) 
where, as in lArnaud et al.l (l2009h . we defined 

Cxsz = ^^ = 1.416 X 10- -Mf^- (B.5) 

The mass weighted temperature is defined as: 

£'""drpg(r)Tg(r)4nr^ 



MW 



(B.6) 



and the factor A = Tmw/Tx takes into account for the difference 
between mass weighted and spectroscopic average temperatures. 
lArnaud et al.l 62009.) find A ~ 0.924. 



